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We consider a charged Brownian gas under the influence of external and non uniform electric, magnetic and me- 
chanical fields, immersed in a non uniform bath temperature. With the collision time as an expansion parameter, we 
study the solution to the associated Kramers equation, including a linear reactive term. To first order we obtain the 
asymptotic (overdamped) regime, governed by transport equations, namely: for the particle density, a Smoluchowski- 
reactive like equation; for the particle's momentum density, a generalized Ohm's like equation; and for the particle's 
energy density, a Maxwell-Cattaneo like equation. Defining a nonequilibrium temperature as the mean kinetic en- 
ergy density, and introducing Boltzmann's entropy density via the one particle distribution function, we present a 
complete thermohydrodynamical picture for a charged Brownian gas. We probe the validity of the local equilibrium 
approximation, Onsager relations, variational principles associated to the entropy production, and apply our results to: 
carrier transport in semiconductors, hot carriers and Brownian motors. Finally, we outline a method to incorporate 
non linear reactive kinetics and a mean field approach to interacting Brownian particles. 
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1: Introduction 

The ubiquitous Brownian motion [T)-[8] remains an outstanding paradigm in modern physics [9]. The theoretical 
framework (TUKU] consisting of Kramers equation (a Fokker-Planck equation in phase space), Smoluchowski equation 
(asymptotic or overdamped contraction of the latter) and the associated stochastic Langevin equation, have been 
widely applied to diverse problems, such as: Brownian motion in potential wells, chemical reactions rate theory, 
nuclear dynamics, stochastic resonance, surface diffusion, general stochastic processes and evolution of nonequilibrium 
systems, in both classical and quantum contexts. More recent applications include thermodynamics of small systems, 
molecular motors, chemical and biological nanostructures, mesoscopic motors power output and efficiency, and the 
fractional Kramers equation applied to anomalous diffusion. Concerning reviews and applications, we mention some 
representative but by no means an exhaustive list of references [10] -[60], some monographs [61] -[70] in addition to the 
"founding papers" [3]- [5J. In his celebrated 1943 Brownian motion paper [IT], Chandrasekhar outlined the method for 
solving a Brownian particle in a general field of force. It took approximately sixty years to report exact solutions for 
the Brownian motion of a charged particle in uniform and static electric and/or magnetic fields |71]-|85] (see also 
some previous related works [551 \S7\ ). 

In Section 2 we generalize and extend previous work |75l 177] and consider a charged Brownian particle in general 
field of force (including magnetic fields), in an inhomogeneous medium (a nonuniform/nonisothermal bath temperature 
profile [8"5]- |113] ). Furthermore we include a reactive linear term (BGK like |114) - |117] ) and obtain a general associated 
Kramers equation. In Section 3 we expand the solution generalizing a previous result |75j yielding recursive relations 
(for other expansions in the literature see for example (64] [73] and |118j - |120] . In Section 4 we use these recursive 
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relations to obtain the Hydrothermodynamical picture of Brownian motion. In Section 5 we present some straight 
forward applications namely, the linear Onsager-like expansions for the entropy production for Brownian motion, 
generalized Shockley carrier transport equations, hot carrier transport and chemically reacting Brownian gases. Finally 
in Section 6 our concluding remarks include further applications of our results. 



2: Generalized Kramers equation for a Brownian charged particle 



Consider a charged Brownian gas composed of N particles (mass m, charge q) in an inhomogeneous bath temperature 
profile, under the influence of external fields (not necessarily uniform in space, nor time independent). Our starting 
point is Kramers equation for the one particle distribution function P(x, v,i), at position x, with v and at time 
t, being in contact with an inhomogeneous bath temperature T(x) (natural units ks = 1) and under a general field 
of force: a mechanical contribution F mcc = —WW, with W — — xF ext (t) + U mec (allowing for potential derived 
forces and non potential homogeneous external forces); an electric contribution F c i oc = qE = —qV(f) and a magnetic 
contribution (Lorentz 's velocity dependent force) F mag =igv x B. The total external force can be casted as a potential 
derived force plus a velocity contribution 

F(x, v, t)= F mec (x,t) + F elec (x,t) + F mag (x,i) - F pot -mw x v (1) 

where 



Fpot = -VU 

U = -xF ext (i) + U mec + q<t> (2) 

w = — B 

mc 

For a charged Brownian particle, Kramers equation for the distribution P(x, v,t) reads (see for example [7511751177] ) 

dP dP 1 . ,dP 18/ T(x)8P\ d ( , „ dP\ 
- + v- + -F x,v,i - = -- V P+ ^i- = — r-VP + r— 3 
at ox m ov t ov \ m ov J ov \ ov J 

where the Brownian collision time is denoted by r, the mobility by A, defined by the relation m\ — r (7 = A -1 is 
Stokes friction coefficient). The mobility and the diffusion coefficient D, are both related to the bath temperature via 
the celebrated fluctuation dissipation theorem D = A/csP (Sutherland-Einstein relation [H[S]). Notice that in our 
scheme the mobility (or r) may also have a position dependent profile; and T is related to fluctuations, as given by 
D = r 2 r. The left hand side of the previous equation may be denoted as the streaming operator L, operating on P 
along the (Newtonian) trajectory v =x, rav = F. The right hand side is denoted as the Fokker Planck collision 
kernel operator Kpp where the first term is the contribution of the (dissipative) Stokes force Fg = —7V and the 
second term the fluctuating contribution, so in explicit form we have 

Kj .,(A,r) = A( T -i v + r £) ( 5) 

and we may cast equation ([3]) in the compact form 

L(F)P(x, v,t) = K fp (t, T)P(x, v,t) (6) 
Let us define a tensorial Stokes force [73 [77] as 

F T s= -A^v-mw x v = -M -1 v (7) 
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where the magneto mobility tensor M is defined (when operating over an arbitrary vector V) as 



,V|tVxw+A(w-V) 



In particular we notice the familiar form |77j for the case B =Bz 

1 tuj 
-™ 1 

1+T 2 UJ 



M(t,u>)= \ 2 I -to 1 | (<)! 



It was shown [7H [77] that by defining such a tensorial Stokes force equation © is equivalent to 



L(f7)P(x, v,t) = K FP (M,r)P(x, v,t) (10) 

thus, the streaming operator includes only de (scalar) potential derived force and the magnetic contribution enters as a 
dissipative non diagonal contribution in the collision kernel (for convenience we have denoted L(f/) =L(F pot = — Vi7)). 
Now as in [751 1115] we add a BGK contribution to the collision kernel denoted by Kbgk ■ This collision operator is 
single relaxation time approximation |114| - [TT7] where the distribution P relaxes (collision time To) to a prescribed 
distribution Pq, the latter instantly equilibrated in the velocity coordinates, to the bath temperature, thus 



Po(x,v,t)=n (x,t)/ (v,T(x)) / (v,m,T(x)) = [^f(x) ) cxp [-^j) 

where fg is a (normalized) Maxwellian distribution, uq is a prescribed density profile. For the Brownian gas we 
define particle density as n(x, t) = J G?vP(x,v,i) where N — J dx.n(x,t). In general n (x,i) needs not to satisfy 
N = J dx.n (x,t) (BGK is not a particle conserving approximation). The BGK kernel is given by 



K BG K(n (x, t), r , fc)P(x, v,t) (P(x, v,i) - P (x, v,t)) (12) 

TO 

and represents the simple chemical reaction A < — > Aq, both compounds with particle density n(x, t) and no(x, t), 
respectively, and with forward (backward) rate /c+(fc_), we have 

dvK BGK P(x, v,t) = -fc+n(x, i) + fc_n (x, t) = - (13) 

where we have defined the rates as k + = k- = Tq 1 . Then, Kramers equation for a charged Brownian particles under 
a general field of force in an inhomogeneous medium, including a linear reactive term, is compacted as 

L(C/)P(x, v,t) - K FP (M, T)P(x, v,t) + K BGK (n (x, t),r )P(x, v,t) = KP(x, v,i) (14) 

This is a generalized kinetic equation in the Fokkcr-Planck (FP) context for a charged Brownian particle. The 
novel aspect in our approach is to include on an equal footing the following aspects: An external magnetic 
field (a velocity dependent force ) rarely considered in the Brownian context; an inhomogeneous bath temperature 
profile T(x) and, in addition to the particle conserving FP collision kernel, we add a non conserving particle colli- 
sion contribution (a generalized BGK mechanism), enabling our approach to incorporate chemical reactions and/or 
generation- recombination processes (for the latter we may consider an inhomogeneous collision time profile tq(x)). 



9n(x, t) 
dt 
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3: Recursive expansion of Kramers equation 



We expand Kramers equation (|14l ) in order to obtain a recursive solution. Some expansion methods include |118] - [T2"0] . 
and we review a method outlined in ref. [75] . The expansion is casted as 



P(x,V,t) = V *„ m ;(x,i)$ nmZ (w) W= \/^ — . V% = ^^ L (15) 



n.m.l— 



1 v , T(x) 



2v T ' 



where $„ m |(w) =4>n(wx)<t > m{w y )4>i(w z ) is the product of (normalized) Hcrmite functions [121j . Following [75] and using 
appropriate recursion relations [121] . equation ([M]) reduces to a difference-differential recursive system of equations 
for the \&'s (v has been integrated out). These recursion relations are better displayed in the compact form, and as 
before |75] we allow for a slowly time varying (compared with r) external temperature field T(x, t). 



where 



with V =-§- and 

9x 



k— n a S n o = tuj ■ A*xAZ n + (A*A+rP)Z„ n = (m, ri 2 , ri 3 ) (16) 

to 



R = §i + 7„ + VA ' + "- ( VA + A T +A ' {« A + A- >¥ + 2? f H 



ni+Tt2+n3yi> 

Zn ~ /III l i5 V 



The asymmetric lowering and raising operators A, A* are defined respectively by 



AiZ a = Z ni _i n2TL3 A\Z n = (1 + WjJ^m-f l,ra 2 ,n3 

^2-^n = ^ni,n 2 -l,n 3 ^2 "^n = ( 1 + "2) ,n 2 + l,n 3 (19) 

A^Z n = Z nij n 2>n3 —i A^Z n = (1 + n3)Z ni) „ 2iTl3 +i 



4: Hydrothermodynamics of Brownian Motion 



We now proceed in a standard fashion [122] - [i"25] and define some relevant moments of the distribution P(x, v,t). 
The usual (lower) moments are: particle density ji(x, t), particle flux vector density J^(x, t), pressure tensor density 
II(x, t), kinetic energy density E(x,t), energy flux vector density Je(x,£) and the total energy flux vector density 
Jq(x, t). Some useful auxiliary quantities are also defined, namely: mass density p(x, t), mass flux density J p (x, t), 
charge flux density J g (x, t), scalar pressure p(x, t) and the stream velocity u(x, t). All of the above are defined by: 

n(x,f) = J dvP(x,v,i), p(x,t) = mn(x,t) (20) 

J M (x,i) = J dvvP(x,v,<) = n(x,i)u(x,<) (21) 
J p {x,t) = mJ M (x,i), J g (x,t) = gJjtf(x,t) (22) 



II(x,i) = m [ dvvvP(x,v,t), p(x,t) = ^ Y^IIij(x,t) 

II u (x,f) = m y dv(v — u)(v — u)P(x, v,i) = II(x,t) — /ju(x,t)u(x,t) 



(23) 
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£?(x,t) = l -m J dv |v| 2 P(x, v,t) = |p(x,t) (24) 
J B (x,t) = \m J dvv|v| 2 P(x,v,t) (25) 
Jg(x,t) = y dvv Qm |v| 2 + f/(x, t)j P(x, v,t) = J B (x,t) + C/(x, t)J M (x,t) 



(26) 



These relevant moments can be readily associated with the expansion coefficients Z n , see equation (|18[) 

n(x, i) = Z 0Q0 , 4>ooo(w) = / (v,m, T(x)) (27) 



^100 

J M (x,t) = ( Z ow 1 (2Si 

^001 



2^200 ■Z'llO ^101 

n(x,t)=n(x,f)T(x)+m| Z uo 2Z 020 ^on ) (2!)) 

■^101 -^011 -^002 



5, 



3-^300 + ^120 + 



1112 



J B (x,t) - -T(x)J M (x,t) + m | Z 210 + 3Z 030 + Z al2 ] (30) 

-^201 + %021 + 3^003 



Then, the recursive relations, equation (|16[) , as outlined in |75) can de compacted into a set of balance equations 
(conservation of mass, momentum and energy, respectively) yielding the hydrodynamical picture. For the particle 
density the recursive relations yields 

Qtl ( Qtl \ 

— + VJ A r = -fe f n + fc_n = — (31) 



at ^ u \m. 

a generalized Smoluchowski equation with a dissipative (reactive) term; for the particle flux we obtain 

+ ( ± + -) 3 M = Jm X c-^VC/-!vn (32) 

if casted in terms of the stream velocity, and using equation (|31[) we obtain the Navier-Stokes equation for Brownian 
flow 

^ + (uV) u+-VII u = 1 (-VC/ + F diss ) , F diss - - fA- 1 + ^) u (33) 

at p m \ tqu ) 

other equivalent versions of equation (|3"2")) in terms of the charge flux density and with F oxt = in equation @ are 



^ + J, = a* (E + P H J g x B) -^V[/ mcc -— VII (34) 
at m m 
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J g = <7* A x (r* , lo) (E--\7U mcc -— VIl) - m A(r*, W )^ (35) 
\ q qn J at 

a generalized Ohm 's law |124j where the effective relaxation time r* , the effective conductivity a* and Hall's coefficient 
are respectively given by 

1 = 1 + 1, CT *( x ,t) = ^ ) i? if ( x ,t) = — (36) 
t* t tq m nqc 

Finally we obtain the energy balance equation 

3E 1 3 

— + V J E = - JmVC/ --(E-E ) Eq = —tiqT (37) 

ot t 2 

Thermodynamics is introduced by defining the kinetic temperature [751 1771 1126) (a generalization of the equipar- 
tition theorem, assuming the Brownian particle to have only translational degrees of freedom ) defined as 



when substituted into ([37| yields 



£(x,t) = ^(x,i)0(x,i) (38) 
3 (90 3 "\ 3 1 

-n— + VJ Q = f -0 + UJ VJ M - - -n o (0 - T) (39) 
We define entropy density 5 and the entropy flux density J5 respectively as [731 1751 177] 

S'(x^) = ~y dvP(x,v,t)ln/«P(x,v,t) 

(40) 

J s (x,i) = - /" dvvP(x,v,t)\nnP(x,v.t) 

where the constant k (lnre = —1 + 3 In is chosen such that under equilibrium conditions S — > S cq (751 1127] (for a 
simple ideal gas, where h is Planck's constant). 

3 

S cq (n cq ,T cq ) = n cq Q + In "^^^ . ™Q( T cq) = (~^T T «i) ( 41 ) 

An entropy balance equation is readily obtained |123j . yielding an expression for the entropy production S5 solely 
in terms of the particle distribution 

+ = *s = J dv(l+l n «P(x,v,*))K(P) (42) 

Finally [73j [75] we define the various thermodynamical potential densities, namely: F (Helmholtz), G (Gibbs), 
[i 1 (intrinsic chemical potential), and fi (total chemical potential), respectively given by [127] 

F(x,t) = E(x,t) - @(x,t)S(x,t) 

G(x,i) = F(x,i) +p(x,i) = rc(x,i)//(x,i) (43) 

/i(x,t) = //(x,t) + f/(x,t) 
the equilibrium expression for the chemical potential is [127] 
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Meq(n eq , T cq ) = -T eq In (44) 

^eq 

All moments, or equivalently, the expansion coefficients Z n , n^O, see equation (1181) . are expandable in powers of t, 
when solving the recurrence relation ()16|) . and considering tlo as an independent parameter. Hereafter, when referring 
to expansion order, we mean in powers of r, for example, the magneto mobility tensor M in equation l[8]). is a first 
order quantity. We further expand the expansion coefficients as 

00 

Z n = Y,T kz( n k) (45) 
fe=l 

For I VT| ^ we have Z [ n k) = for k < integer [^±2] . For the isothermal case |VT| = |£ = we have Z {k) = for 
k < N. 

Now, we discuss several relevant cases concerning the heat bath interaction with the Brownian gas (a gas of real 
particles or some collective property of a given macroscopic system). In order to simplify the discussion we denote the 
Brownian gas as the solute and the heat bath as the solvent (the fluid in which the Brownian particle is immersed). 
As the solute fluctuates and dissipates in the solvent, towards an equilibrium (or non equilibrium steady) state, the 
solute temperature 0(x, t) evolves according to the Brownian dynamics (given by equation [3]) adjusting itself to the 
solvent temperature T(x,i) and the external fields acting on the solute. On the other hand the solvent temperature is 
perturbed by the solute temperature, a very weak interaction indeed, since as a heat bath, T(x, t) relaxes instantly (or 
very fast compared with r) to a prescribed heat bath temperature profile Ti n (x), the latter a stationary function to 
be considered as both an initial and a boundary condition. We realize that VT is not controlled at all in the interior 
of the system but only on its surfaces. In the bulk, the solute performs its Brownian evolution, while the solvent 
temperature relaxes (very fast) to the prescribed profile. This is a unsatisfactory state of affairs since in the bulk, 
the control parameter Tl n (x) evolves as T(x,i) perturbed by the solute, at temperature 0(x,i), which in turn evolves 
following T(x,i), thus a clear cut distinction between controlled external (thermal) fields and the resulting (thermal) 
fluxes cease to exist [125] . 

We analyze three simple cases, with no attempt to model the heat bath at this point (see for example 65, 128 ). 

- Strictly Isothermal Bath, with T(x,t) = T in (x) = T R 

- Strictly Stationary Inhomogeneous Bath with T(x,i) = T; n (x), |VT| 7^ 0, and 

- Weakly Interacting Bath, first order in r, where T(x,t) = Tj n (x) only at the boundaries. 

From the recursion relations (|16l) . for the strictly isothermal case, the solute temperature has a second order 
correction (in r), already obtained by us, as an asymptotic solution to a particular isothermal exact solution [77j . We 
present our result, as in [77 , in terms of the stream velocity u(x, t), in turn expressed as a magnetocovariant derivative 



e(x,t) 



T R + -mu 2 (x,t) 



(46) 



where 



u(x,t) 



■Mcxp - 



n (x, t) exp 



(47) 



n(x,t) r \ T R 

This is a generalization of Shockley's expression for hot carriers effective temperature [1291 1130] . Shockley's result 
coincides with our result at zero magnetic field and homogeneous carrier concentration (provide we identify the hot 
carrier with the solute, at temperature 0, and the solvent with the lattice, at temperature Tr). 

For the weakly interacting bath, to first order in r, we obtain 



0(x,t) = T(x,t) + = T(x,t + I) (48) 

identical to equation (7.2) obtained by 126 , provided we identify the solute with electrons, at temperature and 
the solvent with the lattice, at temperature T. 
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5: Applications 



A)- For the strictly isothermal bath case, we present the second order results, from the recursion relations (fT5f and 
by expanding (|40[) . Our results are equivalent to the asymptotic regime for an exact solution |77j : the local equilibrium 
approximation |131] - [T3"5] is not fulfilled [136 ; in particular the entropy density and the chemical potential are given, 
respectively by the expressions 

S(n(x,t),e(x,t)) = 5 cq (n(x,t),e(x,t))-^-n(x,f)u 2 (x,t) 

(49) 

ju(n(x,f),8(x,t)) = /i cq ()j(x,t),e(x,f)) +{/(x,f) + |mu 2 (x,() 
where the temperature O and the stream velocity u are given respectively by equations (|46p and (|47[) . 



B) - For the strictly stationary inhomogeneous bath case, the local equilibrium approximation holds, in particular 
the entropy density and the chemical potential are given, respectively by the expressions 

S(x,t) = 5 eq (n(x,t),T in (x)) 

(50) 

fi(x,t) = ^ eq (n(x,t),T in (x)) + J7(x,t) 

C) - For the weakly interacting bath, again an as in case A) the local equilibrium approximation does not hold, to 
first order (in r) . The particle flux density is given by the (equivalent) expressions 

JjK-(x,t) = -ra(x,t)MW(x,*) -MV(e(x,t)ra(x,t)) 

(51) 

= -n(x, t)MV(f/(x, t) + 6(x, i)) - D(x,f)Vn(x, t), 
where D = 0(x, i)M is the magneto-diffusion tensor. The total energy flux density is given by 

J Q (x, t)= Qe(x, i) + C/(x, i)J J M (x, i) - |n(x, i)6(x, i)M 3 V0(x, i) (52) 

and where we have incorporated, at no cost, an inhomogeneous collision time profile r(x) as well as a space varying 
magnetic field, with the compact definitions 

M fc V = 0(x) (V+r fe (x)V x w(x)+t2(x)w(x)(w(x).V)) 

= i ; ? ^^2 r,(x)^ir(x) (53) 

l + (Tfe( X )w(x)) 2 fc 

M fc (r,a;) = M(^) M 1 = M fc = l,2,.... 

Notice that in equations (|5T|) and ([52]) according to equation (|48| 9(x, t) can be substituted by T(x, t) since the 
magnetic mobility tensor M, equation ([8]), is already a first order quantity. With the last observation in mind, we 
rewrite the "fluxes" in terms of the "forces", following |125[ 1133] , Thus, we define the matter and heat "forces" as 

Xj*(x,t) = -V (£) X Q (x,i) = vQ^ (54) 
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where 



/i(x, t) = u cq (n(x, t), T(x, i)) + l/(x, i) (55) 
so to first order in r we have the linear flux-force relations 



(56) 



J M (x,t) = LjK«(x,t)X M (x,() +LMQ(x,i)Xg(x,t) 

J Q ( x ,i) = L QM (x,*)Xm (x,i) + Lqq (x,t)X Q (x,i) 
where the coefficients are given by 

L MM (x,t) = n(x,t)T(x,*)M 

L M q(x,£) = Lqa/ (x,i) = n(x, i)T(x,i) (§T(x,t) + f/(x,i) ) M (57) 
L QQ (x,t) - n(x ) t)T( X| t)^T(x l t) + U(x ) f)J M+^(x,i)T 3 (x,t)M 

L+ QM (B)=L MQ (-B) (58) 

the last equation indicating that Onsager relations are satisfied |125j . |131j - [135] . Nevertheless, the Onsager like 
coefficients do not depend solely on the equilibrium values of the state variables (mass density, temperature) . Instead 
they are functions of the external potentials and nonequilibrium state variables, namely n(x,t) and 0(x,i), which in 
turn evolve according to equations (I3I[) and (|39[) respectively (for simplicity of presentation we omit the spatial and 
temporal dependence) 



dn 



—n- 



-VJm — k + n + k-riQ 

(59) 



3 ae /5_„ _ \ . 5_ , 3 1 



^ve + vu ) j M - evj M + (n9M 3 vr) - -— n (e - t) 



2 dt \2 J ni " ' 6 v ' 2t 

and substituting from equation (|48| , within the linear approximation (first order in r) we obtain the coupled equations 
for the evolution of the state variables n(x,£) and T(x,i) (or 0) 



— + — (n - n ) = VMi Y Y =nVU + nVT + TVn (60) 
dt t 

l(" + S)f + 2 rS : = O r + w ) M ' Y + TV(M - Y)+ s v < M »" Tvr > < 61 » 

Equation (jBU)) is a generalized Smoluchowski like advection diffusion reactive equation, and equation (|6ip is a 
Maxwell-Cattaneo like equation, generalizing Fourier's heat equation, incorporates inertial effects [7U1 I126] . |137] - 
[139] . as noticed from the time delay character of equation (f4"8)l . where the time needed for the acceleration of the heat 
flow is considered. 

Thus, even though Onsager relations are satisfied, the coefficient's dependence on nonequilibrium state variables and 
external fields, indicates that any variational principle associated with the entropy production [132, 1341 11351114011141] . 
is not applicable to the present case (69j 1142) 

In the absence of the BGK reactive term, null magnetic field, and for the strictly stationary inhomogeneous bath 
case, we recover the inhomogeneous media advection-diffusion equation, presented by [961 and given by 
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dll{ Q^ = V (A(x) [n(x,t)Vt/(x) + n(x,t)VT(x) + T(x)Vn(x,t)]) (62) 

and for the rigid conductor (constant particle density n), no reactive term and null external fields, we retrieve the 
inertial corrections to Fourier's law, known as the Maxwell-Cattaneo equation |137| - [T3"5] (of the hyperbolic type, 
known also as the telegrapher's equation) 

0T(x,t) , t 9 2 T(x,t) _ K / T(X)t)v 2 T(X)t) + ( V r(x,t)) 2 ) = -VJf(T(x,t)) (63) 
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9* 2 dt 

J F (T) = ^T 3 (x,i)V (^y) = -^(x,t)VT(x,t) = - K (x 1 i)VT(x,i), (64) 

with X = ^A.The right hand side of the last equation is known as the fundamental form of Fourier's law [63] 

D)-Carrier transport in semiconductors, in the Brownian scheme, for the strictly stationary inhomogeneous bath 
case [143, 144 . Consider two Brownian gases, say electrons (density n c , collision time r c , mobility M c = M(t c ,cj c ), 
effective mass m c and charge — e) and holes (density p v , collision time r v , mobility M v = M(t v ,cl> v ), effective mass 
m v and charge e, and with rn c ,v c = w c ,v|eB|).The BGK mechanism is associated in this case to carrier generation- 
recombination with the respective generation-recombination times t® v and with the n°'s the respective equilibrium 
carrier concentrations. Then, by applying equation (|60[) we obtain 



d 1 
—f- = VM C (-en c E+n c VT + TVn c ) - — 



^ - VM v (+ep v V+p v VT + TVp v )--^(p v -p° v ) 



(65) 



Shockley's equations (unidimensional, isothermal, null magnetic field case) are readily recovered [143[ 1144] . Inho- 
mogeneous generalizations of Shockley's equations are of current interest |145j - |152] , 

E)-Our scheme can readily be generalized for several Brownian gases (a = 1, interacting via non reactive 
potential U a /3 and with the usual BGK reactive term associated to each gas. By applying a mean field approximation 
scheme to a many component Brownian gas [153, l54j . as worked out in several contexts in |155] - [16Tj . Kramers system 
of equations can readily be generalized to 

Mt/° ff )P«(x,v,i) = (K FP (M Q ,rj +K| GK )P Q (x,v,i) (66) 

where 

Uffat) = U Q {x,t)+J2 J dyU aP ( X ~y)n p {y,t) n p (y,t) = J dvPp(y,v,t) (67) 



F)- We briefly sketch how to incorporate several chemical reactions schemes, beyond the BGK scheme (deterministic, 
stochastic, photoreactions, etc., see for example [162j - |175] ) in order to model molecular motors. We present two 
chemical reactions that can be approximated at the Kramers equation level, yielding the correct Smoluchowski- 
reactive equation. 

For the nonlinear reaction 

k+ 

A + B ^ C + D (68) 
k. 

the kinetic equations are retrieved if we start consider the BGK like scheme 
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~k + n B P A + k^ncn D fo(v,m A ,T) 



-k + n A P B + k-n c n D fo(v,m B ,T) 



-k^n D P c + k + n A n B f (v,m c ,T) 



-k-n c PD + k + n A n B f (v,m D ,T) 



(69) 



Finally, for the reaction cycle 



we reproduce the correct kinetics with the approximation 

K BGK P a(*,v,t) = -k a P a + k a -ifo \v,m a ,T)n a _i 



n = m 



a = I,. J 



(70) 
(71) 



6: Concluding remarks 



We have extended our previous work on charged Brownian particles |73l 1751 177] , in order to obtain a consistent 
expansion scheme in powers of the collision time. We presented the complete hydrothermodynamical picture for 
charged Brownian particles in the Kramers equation scheme, considering the action of external magnetic, electric and 
mechanical fields, chemical transformations within the BGK scheme, and furthermore space dependent thermal fields 
(inhomogeneous media). 

We developed a recursive method, in order to expand in powers of the collision time r, the several moments of the 
solution of Kramers equation, enabling us to compute the governing equations for mass (charge) , momentum (stream 
velocity) and energy (heat) flow. In the appropriate limits we retrieve previous results, and present novel formulations 
and results with immediate physical consequences. 

Applications A, B and C (previous section) clearly indicates that the Brownian scheme cannot be modelled, in general 
and in earnest, by the local equilibrium approximation, no variational principle associated with entropy production 
is applicable, and hyperbolic extensions of Fourier heat conduction law must be considered even in lowest order. 
Therefore, variational principles relevant to the Brownian scheme, must go beyond entropy production considerations, 
and include among others, free energies, efficiency factors and power output as well |691 I70j and connected to 
nonequilibrium fluctuation and work relations |176j - [l85] . 

In the following Applications items we briefly outline, from the formalism developed in this work, some novel schemes: 
D) the incorporation of magnetic field effects on the nonequilibrium state variables (including an application for carrier 
transport in inhomogeneous semiconductors) ; E) inclusion of interacting Brownian fluids via mean field approximation 
schemes, and F) inclusion of chemical reactions to Kramer's approach to Brownian motion. 

Our future work concentrates in variational principles associated to the Brownian motion, inclusion of chemical 
reactions and the magnetic field effects in hot carrier transport in semiconductors |186j . 
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